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Abstract. - Last year, Blanco and Fournier (Blanco S. and Fournier R., Europhys. Lett. 
2003) calculated the mean first exit time of a domain of a particle undergoing a randomly 
reoriented ballistic motion which starts from the boundary. They showed that it is simply 
related to the ratio of the volume's domain over its surface. This work was extended by 
Mazzolo (Mazzolo A., Europhys. Lett. 2004) who studied the case of trajectories which start 
inside the volume. In this letter, we propose an alternative formulation of the problem which 
allows us to calculate not only the mean exit time, but also the mean residence time inside a 
sub-domain. The cases of any combinations of reflecting and absorbing boundary conditions 
are considered. Lastly, we generalize our results for a wide class of stochastic motions. 



Introduction. - Recently, motivated by the study of diffusive trajectories performed by 
animals in various conditions, Blanco and Fournier [1] reported an important result concerning 
the mean first exit time of Pearson random walks [2] in a bounded domain. Beyond the 
applications for animal trajectories, numerous physical systems are concerned, such as the 
neutron scattering processes [3,4]. Pearson random walks can be defined by the trajectory of 
a particle which is submitted at stochastic times to random reorientations of the direction of 
its constant velocity v. The stochastic times are exponentially distributed with mean value 
1/A. For this category of random walks, Blanco and Fournier showed that the mean first exit 
time < < > of a random walk starting from the boundary of a finite domain is independent of 
the frequency A and is simply related to the ratio of the domain's volume V over the surface 
S of the domain's boundary, namely in three dimension: 

4 V 

<t>^-^ (1) 
V b 

The fact that this result does not depend on the frequency A but only on the geometry 
of the system seems counter-intuitive, although Blanco and Fournier propose some heuristic 
explanation for it. Very recently, following the demonstration of Blanco and Fournier, Mazzolo 
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[5] extended their results and gave a relation between the nth moment of the first exit time 
of a particle starting inside the domain, and the n + 1th moment of the first exit time of a 
particle starting from the boundary. 

In this paper, we present an alternative approach, which allows us to extend these results 
in three directions. First, we consider the case of any combination of absorbing and reflecting 
boundary conditions, which is useful in many applied stochastic processes. Second, we study 
not only first exit time properties, but also residence time within a sub-domain properties 
(prior to the first exit of the domain), which is a more general quantity. Last, we show that 
these results can be generalized for a wide class of stochastic motions. 

First exit time. - Before considering the case of a general motion, we restrict our discus- 
sion to the previously examined Pearson walks, which corresponds to a succession of determin- 
istic ballistic movements in the bounded domain V delimited by an absorbing boundary S, 
interrupted by instantaneous and isotropic redistributions of the velocity v. Let p{r' , v' , i|r, v) 
be the conditional density probability at time t that the particle is at position r' with a ve- 
locity w', given that it starts initially at position r with a velocity v. This quantity satisfies 
the well-known backward Chapman-Kolmogorov differential equation [6] : 



dtp{r' ,v' ,t\r.,v) = V ,pp{r' .,v' ,t\r,v) ^ / dv''[p{r' ,v' ,t\f.,v'') — p{r' ,v' ,t\r,v)\, (2) 



d 

where ad = ■ildT stands for the solid angle in c?-dimension and the integral f dv is taken over 

-'-12^ 

unit vectors in the direction of the velocity. This equation can be easily converted into an 
equation for the moments of the first exit time T(r, v) starting from r with the initial velocity 
v[I\. 

v-Vptn{r,v) + — [ dv"{tn{r,v") -tn{r,v)) = -ni„_i(f, t/), (3) 

where tn{r,v) is the nth moment of the first exit time, which verifies tn{r,v) = on the 
absorbing boundary (r S S), for all n > 1 and for an outward velocity. Noting that the 
symmetric quantity J dvdv"tn{r,v") ~ J dvdv"tn{r,v) obviously equals zero, and that v ■ 
^pt„(r, i/) — div[tn{r,v)v), the integration of equation (|3J) over all possible initial positions 
and velocities gives: 

d-D / drdiv{tn{r,v)v) — ~n dv dftn~i{r,v). (4) 



Applying the Gauss divergence theorem on the left-hand side of the equation Q), we can write 
that: 

<tn>S^Vd-^ <tn-l>V, (5) 
V Zj 

where the two averages involved in equation ((SJl are defined for any function / as follows : 

< / >s= ^ [ dv [ dS • vf{f,v), and < / >y= -J— ( di ( dff{f,v) (6) 

vT,ad J Vad J Jy 

d-l 

2tt— 

with ad — , -I is the inward flux of a unit, isotropically distributed vector through 

(d-i)r(i^) ^ y J 6 

Cd r 

a unit surface, and 7]d = — = ^/ntd—l) . ^. is a dimension dependent constant. Equation 

ad r(|) 
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Fig. 1 Fig. 2 

Fig. 1 - A particular particle's trajectory which starts and ends on the absorbing boundary Sg^j^g, for 
an absorbing region inside a closed reflecting surface E 

Fig. 2 - A particular particle's trajectory for a reflecting box with a hole (aperture) Sg^j^g 



(jSJ can be rewritten as: 

< h >s= — and \/n > 1, < <„-i >v= : (7) 

which corresponds to the results previously obtained by Blanco and Fournier [1] and Mazzolo 
[5]. Parenthetically, we know that since < ^2 >e>< h >|, iQ yields a general lower bound 
for < ti >v- 

<h>v>-^<h>.^l^. (8) 

General boundary conditions. - We now extend these results, which take unexpected 
simple forms, to the general case of mixed boundary conditions. The domain's boundary 
E is now supposed to be composed of any combination of absorbing and reflecting parts 
(see fig lll2() . Let Eg^j^g be the absorbing part of S, and ^^.g^ its reflecting part, so that 
E = Eg^]-|g + ^refl' assume that initially the particle starts from the absorbing surface. 
Whether we consider that the reflecting boundary is perfect, i.e. that the inward velocity 
angle equals the outward velocity angle, or whether we consider that the particle's velocity 
angle is uniformly redistributed after each collision with the reflecting boundary, we have 



[dv [ 



dE-z7i„(f,z/) = 0. (9) 

refl 

In this case of mixed boundary conditions, equation lO becomes : 

<ti>s, --^^ , and Vn > 1, < t„_i >v= ^ (10) 

Let us, for instance, consider animals leaving their home through its boundary E^^j^g, and 
exploring a domain V (see Fig^. The previous equation gives a simple estimate of their mean 
return time to home if the boundary Ej.gg is reflecting for them, provided that the interrupted 
ballistic motion correctly models their behaviour. 
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Splitting probabilities. - Consider now a closed volume V with two absorbing surfaces 
El and S2, and let us compute the absorption probability Ili{f,v) on the surface Si for a 
particle starting initially from the boundary Sg^^^g = Ei + S2. This probability obeys the 
differential equation [8]: 

V- ni(f, v) + — [ d?7"[ni(f, w") - ni(f, v)] = o, (u) 

CTd J 

for which the boundary conditions, for an outward velocity are ni(r, iT)^?^^! = 1 and 
Ii-i{r,v)pe-E2 — 0- Following the same lines as previously, one can show that 

< >-abs- E^- 

In other words, the splitting probability through absorbing boundary is uniform, when starting 
uniformly from any point of it. 

Residence time. - A relevant quantity involved in the theory of chemical reactivity and 
in other similar problems [9-12], is the residence time in a sub-domain V' (see figH)) of the 
domain V. In order to study it, we assume that the particle can disappear in V' (but not 
outside V) with a constant and uniform reaction rate k. The central quantity is now the 
survival probability of the particle with respect to the reaction in volume V (i.e. considering 
that an absorption on the boundary is not a killing process), which can be expressed in two 
ways [12-14]: 

S{t\r,v) = [ dr'di'p{r',v',t\r,v) =< e-'''^''-''^'^ > (13) 



where Tt is the residence time up to the observation time t inside V' of the particle, starting 
from {r,v), and p(r', z;', t|r, ■u) is the conditional density probability of the particle which 
satisfies the following backward Chapman-Kolmogorov differential equation: 

dtp{r' ,v' ,t\f,v) — V pp{r' ,v' jtlr'jv) -\ f dv"[p{r' ,v' ,t\f,v") — p{r' ,v' ,t\f,v)] 

J 

-klv'{r)pir',v',t\f,v) (14) 

where ly'ir) is the indicator function of V' which equals zero if r ^ V and equals one if 
r S v. Note that the boundary conditions give S{t\f,v) = 1 for an outward velocity of a 
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particle on the absorbing surface Sg^|-,g. Indeed, such a particle absorbed on this surface will 
never disappear within the sub-domain V' . Writing down the equation H14|) for the conditional 
survival probability, we obtain: 

dtS{t\f,v) = v.^pS{t\f,v) + — [ dv"[S{t\f,v^')- S{t\f,v)] 

CTd J 

-klvir)S{t\f,v). (15) 

As S{t\f, v) is a bounded non-increasing function of t, it tends to a finite limit as < — > oo, so 
dS 

much — > 0. Letting t ^ oo in 1151) . integrating it over r and u, taking into account 

at 

the boundary conditions and using the expansion (see H13|l when t — > oo) 

5oo(f,t^) = 5^(-l)"^T„(f,z;) (16) 

n=0 

where r„ (r , v) is the nth moment of the residence time of the particle inside V' for an infinite 
observation time, we derive our main result for the residence time: 

< Ti >s , = ^- and Vn > 1, < r„_i >y,= (17) 

« Sabs n<n >E^^^g 

Note that the mean fraction of time spent inside V' prior to the first exit of V is given by 
< Ti >Sa^|^g / <ti '*s^|_^g= V /V, which can be seen as an ergodic type property [8]. 

Generalisation. - Let us now assume that the particle undergoes a more general stochas- 
tic diffusion process, interrupted at stochastic times by an instantaneous redistribution of its 
velocity, without changing its position. The conditional density probability p{r' ,v' ,t\f,v) to 
be at position r' with velocity v' at time i, starting from (r, v) at time 0, corresponding to 
the general diffusion process obeys the backward equation 

F 

dtp = VidxiP + —dy^p + Dijdl p = L+p (18) 
m ' 

where m is the mass of the particle, F{r, v) is the force exerted on the particle in state 
(r, ■y), Dij is the general diffusion matrix, and L"*" stands for the adjoint of L. In many 
situations of physical interest, F is the sum of a conservative force due to a potential J7, 
—dpU, and of a friction force —r]v due to a surrounding fluid, rj being a constant friction 
coefficient. Furthermore, Dy = SijDrj^ /m? ^ D being the usual diffusion coefficient of the Fick 
law, which is independent of r in an homogeneous medium. Now, we assume that this general 
diffusion process is interrupted at stochastic times, for which the velocity v of the particle is 
instantaneously redistributed with a given transition rate q{f,v'\f,v). The overall stochastic 
process, including the velocity reorientation, obeys the equation: 



dtp{r' ,v' ,t\r,v) — L+p(r', u', i|r, ?j) -I- y dt;"(7(r, u"|r, ?j)[(p(r', w', i|r, ti") — p(r', tj', i|r, z/))] 

= L'^p{r' ,v' ,t\r,v). (19) 

The conditional probability density at time p(r', w', i|r, v) also obeys the following forward 
equation: 

dtp{r' , v', t\f, v) — Cp{r', v' , t\r, v). (20) 
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Let US assume that this equation, when the whole boundary E is assumed to be reflecting, 
admits a stationary solution po{f, v). For instance poif^ v) is the thermodynamic equilibrium 
^-{u+^^)/kBT j£ equation (|18|l is the usual diffusion equation of a particle in an equilib- 
rium fluid at temperature T and if the velocity redistribution process defined previously does 
not affect this thermodynamic equilibrium, which is the case if ^ dv' {q{r,v'\f,v)pQ{f^v) — 
q{r,v\f,v')pQ{f,v')) — 0. Then it is possible to generalize the previous results. Indeed, let 
us consider the residence time of the particle in a sub-domain V' G V , and let S^^j^g be the 
absorbing part of the total surface S. As previously, we assume that inside V , the particle 
can disappear with a constant and uniform reaction rate k. Its conditional survival probabil- 
ity at time i, S{t\r,v) — J dr' dv' p{r' , v' , t\r, v) , starting from (r S V,v) obeys the backward 
equation 

§iS{t\r,v) = C+S{t\f,v) - klv'{r)S{t\f,v) (21) 

where lv'{r) = 1 if r G V and otherwise. Following the method used for the Pearson 
random walks, we let t — > -|-oo in (|21|l . and we multiply this equation by the equilibrium 
probability po{r,v) computed for the reflecting conditions on the boundary S of V. The 
integration of equation (|21|l over all r S V and v, leads to: 

/ ,f[,,Mf.^C*S^(f.v)^k[ df/d*.K.7)5„(f,,l. (22) 

JrGV J JreV J 

The left hand side of (|22(l can be transformed into 



dr dvSoo{r,v)Cpa{f,v)+ / dS • dvvSoc{r,v)po{r,v). (23) 
f£V J JfeT. J 

The first integral vanishes, since Cpo{r, v) = if r (z V . In the second integral we separate the 
inward and outward integration over v in the surface integral. The inward part can be written, 

using the expansion itTBI) of Sqo'- > (— 1)"'^^— T'^ < Tn >e , , where J — — f-^v dE • 

^ n! abs '^'^^^abs 

n— 

/ dvvpQ{r,v) is the one-way equilibrium probability current on Sg^|-,g, and the surface average 
of a quantity f{f,v) is defined as: 

<f>^.h.= --Tf / dvvf{f,v)po{r,v). (24) 

""^^ -^^^^^^abs 

The outward part is actually equal to J, as Soo{'r,v) = 1 if r G S and v points outward. 
Then, comparing both sides of (|22|l . we generalize the relations (117(1 for the moments of the 
residence time within V': 

<^i>^abs=7 ^''^ <^">^abs="^^^^T^ 

where P — Jp^y, drdvpo(r, v) is the equilibrium probability of volume V' . 

Considering now the special case V' = V, we obtain the relations for the moments of the 
first exit time from V: 

1 1 < >s , 

< ti >s , = ^, and < tn-i >v= : —■ (26) 

abs J 71 < <i >s , 

abs 
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which generahze the formulas Q . In the same way, the resuhs of the sphtting probabihties 
H12|l can be generahzed for the general stochastic motion discussed here. Equation H2t)|) takes 
a particularly simple form if we now consider a fluid of identical particles and if we replace 
Po{r, v) by the equilibrium density of particles in the phase space (r, v). Then, the probability 
current J is replaced by the equilibrium one-way flux of particles (j) through S, and equation 
becomes: 

< ti >s= (27) 

A'^ being the number of particle in Y . This formula can easily be understood intuitively: if at 
time the boundary is supposed to be reflecting, the number of particles leaving V per unit 
time should be ^ <:ti>s ' that the equilibrium is maintained inside if y is supplied with 
an entering flux of new particles (j) ^ <t^^- 

Conclusion. - We have shown that the geometrical relations previously obtained for Per- 
son random walks in bounded domains, are particular cases of very general relations between 
residence times for a large class of stochastic processes. Thus, they can be very useful in many 
applications, when the evolution equations can not be solved exactly. Furthermore they can 
be extended to intermittent systems, which are frequent in biology [15]. Lastly, with relevant 
changes, this work can be extended to discrete space and time systems. Such extensions are 
in progress. 

* * * 

We thank Sylvain Condamin for useful discussions. 
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